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1. Introduction 



Many fundamental models of statistical mechanics can be expressed as 1 + 1- 
dimensional field theories with a boundary. A number of these field theories can be exactly 
solved by applying the powerful techniques of integrability [jl[. Such models include the 
Kondo problem, a number of models of dissipative quantum mechanics, and the model we 
treat here: a free boson with a periodic potential on the boundary. 

We consider the "boundary sine-Gordon" problem of a massless scalar field on 
the segment < a < L. The Lagrangian is 



This model arises in a number of contexts because one can often map interacting fermions 
(the Luttinger liquid or massless Thirring model) onto a free boson. Coupling the fermions 
to a boundary potential, for example, gives this boundary interaction for the bosons f^. 
This model has also appeared in dissipative quantum mechanics, where it describes a 
particle moving in a periodic potential with a frictional force [Q. 

The model in the bulk (a free massless scalar field) is, of course, a conformal field 
theory. For g=0 and g oo, the combination of left and right conformal symmetries 
is preserved in the presence of a boundary However for (3'^ < Stt, the interaction 
breaks the scale invariance, with the parameter g providing a mass scale. At g=0, the 
boundary dimension of the operator cos[/3(^/2] is x = For 0< x <1 a perturbation by 
this operator is relevant and induces a renormalization group flow between two different 
conformal boundary conditions. Our purpose is to study this flow explicitly. 

To understand the boundary condition induced by the interaction, we examine 
the energy-momentum tensor T^jy = ^df^(f)d,y(j), which satisfies the current conservation 
d^T^i, = 0. The Hamiltonian is defined as if = Tttda — g cos[P (f)(0)/ 2]. The inter- 
action adds a boundary contribution to H, so that the boundary itself can hold energy. 
Without the boundary interaction, we would require that the component T^t vanish at 
a = and a = L to ensure that no energy crosses the boundaries (i.e. dH/dt = 0). 
Thus in general, energy conservation gives a dynamical boundary condition, requiring that 
Tat = dtg cos[l3(f}/2] at cr = 0, or 
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This is a classical equation, but the same relation holds in the quantum field theory after 
a renormalization of parameters . Thus the boundary interaction causes a fiow between 
Neumann boundary conditions {d(j){0) = 0) in the UV and Dirichlet (^(0) = Ann/ (3, n 
integer) in the IR. These two boundary conditions are dual to each other. 

To solve this model, we shall make use of the powerful constraints of integrability. 
There are an infinite number of conserved kinematic quantities for the free boson. A generic 
boundary interaction will destroy them. However, a boundary interaction of the form ( |1 . IJ ) 
(or equivalently, a boundary condition of the form (|1.2| ) ) still preserves an infinite number 
of appropriately-modified conserved quantites. This can be seen explicitly in perturbation 
theory and in the classical limit, by taking the massless limit of the sine-Gordon model 
with a boundary f^. 

A particularly useful approach to integrable theories with and without boundaries 
is to find the exact S matrix of the physical quasiparticles of the theory From 
the S matrix one can calculate a number of exact thermodynamic quantities such as the 
free energy, specific heat, susceptibility and boundary entropy. This is done by using 
the thermodynamic Bethe ansatz We describe this technique in the presence of a 

boundary in sect. 4. 

To understand the states of our theory, we look at the system far away from the 
boundary. Here the "bulk" theory is just that of a free massless scalar field. Ordinarily, 
with a massless bulk theory, one describes the excitations in terms of plane waves (perhaps 
around some non-trivial topological background). However, this is not a convenient basis 
to work with once one has included the boundary interaction. The reason for this is 
simple. The plane waves are no longer eigenstates of the conserved charges, so a plane 
wave scattering off the boundary could result in a possibly-horrible combination of other 
plane waves. 

A basis where the states have well-defined eigenvalues of the conserved charges is given 



by the massless limit of sine-Gordon solitons [^. In other words, we consider therefore the 
problem given by the Lagrangian (|1 . 1| ) as the G — ^ limit of a problem which is massive 
in the bulk, with perturbation dC = —G da cos (3(f). With this particular combination of 
bulk and boundary perturbations (notice the (3 in the bulk and the (3/2 on the boundary) 
the model is integrable 0. In sect. 2, we will make this change-of-basis explicit by showing 
that the massless solitons obey the classical equation ( |1 . 2| ) , and that they retain their form 
after scattering. 
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When one changes basis from the free plane waves to the massless sine-Gordon sohtons, 
the scattering of the bulk particles is no longer trivial or even diagonal. We describe such 
processes by massless scattering [|T2| , p!3| , |TT| . The integrability of the model is crucial to 



such a description. In a basis with well-defined charges, all collisions must be completely 
elastic, so that momenta are conserved individually and particle production is impossible. 
Even with this crucial simplification, it is still not immediately obvious what it means to 
scatter two particles traveling in the same direction at the speed of light. There are several 
ways to define a bulk S matrix in the integrable models we are studying. One is as the 
phase which results when commuting two particle-creation operators (the Zamolodchikov- 
Faddeev algebra). This phase is equivalently realized as a matching condition on a two- 
particle wavefunction 'i/'(cri,a2) between the regions o"i ^ o"2 and ai ^ a2. One can 
explicitly compute this phase in a lattice Bethe ansatz by forming a two-particle state and 
bringing one of the two particles around the periodic world, i.e. sending ai ^ ai + L, 
where space is a circle of circumference L . In other words, one studies the monodromy 
properties of massless particles [ p!5| . 

We emphasize that by "massless" limit, we mean that the quasiparticles are massless 
(they have the dispersion relation E = ±P). The theory is not scale-invariant because of 
the boundary interaction. With this scale (we call it Tg), boundary scattering can depend 
on the momentum of the incident particle, because there exists a dimensionless parameter 
P/Tb- Similarly, thermodynamic quantities like the free energy depend on the ratio T/Tb- 

The outline of this paper is as follows. In sect. 2, we discuss the classical limit of the 
problem in order to give a more precise understanding of massless scattering. While it is 
conceptually important, the details are not required for what follows. In sect. 3, we give 
the exact quantum S matrix by taking the massless limit of the boundary S matrix for 
the sine-Gordon model 0. We show that it has all of the desired properties, and that it 
agrees with a previously-derived result at 0^ = An [|16|. In sect. 4, we use the S matrix to 
derive the boundary free energy, thus showing how the boundary entropy flows between 
the Neumann and Dirichlet fixed points. In sect. 5 we present some thoughts on future 
directions. 

2. Understanding massless scattering via the classical limit 

Studying the scattering of massless particles has given a wealth of useful information 
generally unavailable by other means; for example, as we discuss in sect. 4, it enables the 
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calculation of the exact free energy (c-function) flowing between two flxed points (for a 
review, see One subtlety is that in such a description there can be non-trivial bulk 

scattering even if a massless theory is linear and obeys superposition, and thus appears to 
have trivial monodromy. For example, in the Kondo problem, the bulk model consists of 
free massless electrons; in this paper we discuss a free massless scalar fleld. The reason 
for the non-trivial scattering, as indicated in the introduction, is that when coupling to a 
boundary, it is much more useful to make a (non-linear) change of basis of particles. The 
advantage is that boundary scattering is much easier to describe; the added complication 
is that the scattering in the bulk is no longer trivial. Thus in the multi-channel Kondo 
problem, the appropriate quasiparticles are kinks in a potential with mutiple degenerate 
wells whereas in this paper they are kinks (and breathers) in a potential with an 
infinite number of wells. This change-of-basis is often not very explicit, since it can be 
deduced indirectly from the Bethe ansatz solution, or inferred from demanding the appro- 
priate symmetry structure. Our purpose in this section is to elucidate massless scattering 
and hopefully give some physical insight by describing the classical limit of an important 
example. 

We consider here a classical scalar field (^(a, t) satisfying the Klein-Gordon equation: 

- dl<P = 0. (2.1) 

A dimensional parameter g is introduced by requiring that theory live on the half-line 
[0, oo), with boundary conditions: 

M.=o = 9sm{^i<p-cPo))l^, , (2.2) 
for some constants g and (po- (Classically, one can scale out the parameter /3 in (|1.2| ).) 

2.1. The massless limit of sine- Gordon 

In many massless scattering problems an intrinsic bulk massive perturbation has been 
implicitly or explicitly introduced, and then set to zero. Moreover, the perturbation is 
always chosen so that the theory remains integrable off the critical point. The fact that one 
is considering a massless limit of a massive integrable model leads one to a preferred basis 
of 'wave packets' or 'wavelets' for the massless theory. This enables one to give a particle 
interpretation to the massless theory, and defines the concept of massless scattering. A 
given massless theory can, of course, be the limit of several massive integrable theories 
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(e.g. the Ising model with a magnetic or energy perturbation), and as a result there can 
be several different massless scattering matrices, leading to distinct descriptions of the 
massless theory (see, for example, [jl^). A bulk mass can be introduced to our example 
by considering ( |2.1|) to be the m — + limit of the sine- Gordon equation: 

d'^cj) - diet) = -m^ sm{(f)) . (2.3) 

This massive model was shown to be integrable in the presence of the boundary condition 
(P) in i- 

For either of these partial differential equations we want to consider wave packets 
of finite energy (which means that (p must approach a constant value suitably fast as 
a —>■ ±00). For any non-linear, integrable, partial differential equation, there are infinitely 
many conserved quantities, and in the limit in which the mass vanishes, these conserved 
quantities become conserved quantities for the massless theory. For the sine-Gordon theory, 
the massless limit of these conserved quantities are special polynomials in the partial 
derivatives of (p. The most general characterization of the preferred wave packets in the 
massless theory are those wave packets that are eigenfunctions of the massless limit of 
these conserved quantities. An equivalent characterization of these massless wave packets 
is to take the non-dispersive wave packets of the massive theories and take an appropriate 
limit in which their energy remains finite when m ^ 0. 

There are two types of finite-energy solutions of the classical sine- Gordon equation 
( p.3|) : solitons, which are time-independent and topologically non-trivial, and breathers, 
which are time-dependent and topologically trivial. Intuitively, a breather can be thought 
of as a bound state of a kink and an antikink oscillating in and out (i.e. breathing). In this 
section, we will discuss only the solitons; the analysis for the breathers follows analogously. 

A major triumph of the theory of non-linear partial differential equations was the con- 
struction of explicit solutions of ( p.3|) for any number of moving solitons (see, for example, 
||19|| ). The solitons' energies and momenta are conveniently expressed in terms of rapidities 
ctj, defined by Ej = mcoshctj and Pj = msinhctj. The velocity of each is thus given by 
tanhctj (positive for a right-moving soliton). We have set the speed of "light" to be 1. 

To find the scattering of two solitons in the massless limit, consider a two-soliton 
solution of (|2.3|) on (— cxd, 00). This solution is usually expressed as: 



((7,t) = 4 arg{r) = 4 arctan ( ) , (2.4) 



where the r- function solution is given by: 



r = l-eiesftanh^^^^ E,ia-a)-E,ia-b)+p,t+p,t 

V 2 ; (2.5) 

The constants a and b represent the initial positions of the two solitons, and ej = +1 if 
the j^^ soliton is a kink, while ej = —1 if it is an anti-kink. 

For a wavepacket to have finite energy in the massless limit m — 0, the rapidity |a| 
must go to infinity. We thus define a = A + 9, and let A ^ oo such that the parameter 
/i = ime^ remains finite. The energy and momentum of a right-moving "massless" soliton 
then reads E = P = /ie^. For a left mover, a = —A + 9, and its energy and momentum 
read E = —P = iie~^ . 

Suppose that both of these solitons are right-moving. Then the massless limit yields: 

r ^ I - ,^,^^-^^-E.iv-a)-E,ir,-b) ^ ^| eie'^^^''-") + e2e-^-^^-'^ ] (2.6) 
where rj = {a — t) and 

A = - log[(tanh(0i-^2))']- 

This is manifestly a solution of the Klein-Gordon equation, and by construction is an 
eigenfunction of all of the conserved quantities of the massless limit of the sine-Gordon 



equation. (It is shown in [|TT[ that it remains an eigenfunction after quantization — the 
expectation values {E'^'^~^) are precisely the quantum conserved charges derived in [pO| .) 
It is almost a superposition of two single-soliton wave packets. Observe that: 

arg\l - 6163 e--^i(''-")--^2(^-b) ^ ^| ^^^-E,{-n-a) ^ ^^^-E^{n-b) |i 

(2.7) 

The factor A thus measures the extent to which the two-soliton solution is not a superpo- 
sition of one-soliton solutions. 

More precisely, consider the limit a 00, rj — »• 00 so that Ei{r] — a) is finite. This 
corresponds to moving the first kink off to cr = +00 and following it. The r function 
collapses to the one-kink form t = 1 + ieie~^^^'^~°'\ Moving this soliton through the 
second one corresponds to taking it to a = —00, or taking the limit a —00,77 —00 
(with El [rj — a) finite) . Discarding an overall multiplicative factor (which is irrelevant in 
the computation of ^ = 4arg(r)), we see that in this limit, r = 1 + ieie~^^^'^~°'^~^ . Thus 
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these preferred Klein-Gordon wave packets exhibit non-trivial monodromy. The foregoing 
time delay A is precisely the classical form of a massless scattering matrix. One can then 



semi- classically quantize this [ pT|p2| ; the analysis is identical in massive and massless cases, 
because the difference 6i — 62 is the same in both cases. We will bypass this step because 
the exact quantum answer is already known [Q. 

One obtains the same A for two left-moving solitons. For a left-moving and a right- 
moving soliton colliding one easily sees that the massless limit of (tanhi(ai — ^2)) is 
unity. The solution collapses to the superposition of a left-moving wave packet and a 
right-moving wave packet exactly as in ( p.7|) , with no time delay. This is the classical 
manifestation of the fact that the left-right quantum scattering matrix Slr elements are 
at most rapidity-independent phase shifts. 

2.2. Klein-Gordon with a non-linear integrable boundary condition 

Consider now the Klein-Gordon equation on [0, 00) with the boundary condition ( |2.2| ). 
The non-linear boundary conditions destroy the naive conserved quantities of the Klein- 
Gordon equation. However, by reversing the argument in the appendix of one finds that 
the system will still have higher-spin conserved quantities, namely those of the massive sine- 
Gordon system. Thus the bulk theory is Klein-Gordon, but in this instance there are a 
preferred set of wave packets that scatter from the boundary without dispersion. These 
wave packets are once again the massless limit of the sine-Gordon solitons. 

A more direct way of seeing the integrability of the Klein-Gordon (and indeed, sine- 



Gordon) equation with boundary conditions ( p.2|) is provided by the analysis of [|3l. The 
idea is to show that the method of images can be used on (—00, 00) even in the non- 
linear system, so as to replicate the boundary conditions (p.2| ) on [0, 00). The scattering 
of a kink, or anti-kink, from the boundary can be described by a three-soliton solution 
on (—00,00). These three solitons consist of the incoming soliton, its mirror image with 
equal but opposite velocity, and a stationary soliton at the origin (to adjust the boundary 
conditions). If one takes the infinite rapidity limit of this three-soliton solution then the 
stationary soliton simply collapses to an overall shift of (/> by a constant, while the mirror 
images (since they are moving in opposite directions) reduce to a superposition of two 
wave packets. One thus obtains: 

(f, = (Po + 4 ar^[l + zeie-^(^-")] + 4 arg[l + ie2e-^^'^-^^] , (2.8) 
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where ^ = a + t, r] = a — t. By direct computation one finds that this solution satifies ( |2.2| ) 
with: 

The constant = —E{a + b) represents the delay of the reflected pulse. If one defines 
the classical boundary scale 9b via g = 2(Ue~^^, then this delay may be written as 

Ab = log(-eie2 tanhi(^-0B)) • (2.10) 

Note that the sign €2 is to be chosen so as to make the argument of the logarithm real. 
This determines whether the reflection of a kink will be a kink or an anti-kink. Thus we 
see that 9b is the scale at which behavior crosses over from the region of the Neumann 
critical point (where the classical boundary scattering is completely off-diagonal) to the 
Dirichlet boundary critical point (where classical boundary scattering is diagonal). This 
result is completely consistent with the infinite-rapidity limit of the delay computed in 

Before concluding this section, we wish to observe that the foregoing is not simply a 
property of some solitonic solutions, but is a general feature of the sine-Gordon equation 
and its massless limit. The sine-Gordon equation (and any other integrable equation) 
admits Backlund transformations. In particular, consider 

d,i<P + i;) = mD sin(^) ; d^{cP " ^) = ^ (^) ' ^^'^^^ 

where D is an arbitrary parameter. The point is that (p satisfies the sine-Gordon equation 
if and only if i/j does so as well. Taking m while keeping mD finite, the Backlund 
transformation becomes: 

a^(</, + tp) = mD sin (^y^) ; d^i(p - tp) = . (2.12) 

One can easily check that this transformation implies that (p satisfies the Klein-Gordon 
equation if and only if tp does so as well. However, suppose that tp is a solution to Klein- 
Gordon on [0, 00) satisfying Dirichlet boundary conditions cp = cpQ. It follows immediately 
from (|2.12| ) that tp is a. solution to Klein-Gordon satisfying ( |2.2| ) with g = ~mD. For any 



left moving wave-packet, /(^), the Dirichlet scattering solution (p is given hj (p = (po+f{$,)~ 
f{—r]). The scattering solution, tp, for the same incoming wave-packet reflecting from the 
boundary conditions (|2.2| ) can then, in principle, be found by by solving ( |2.12|) . That is, 
if one writes tp = /(^) + r{r]) then r satisfies r'^r]) = f'{r]) + g sin((r(?7) + f{—f]) — (po)/2). 
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3. The quantum S matrix 

We find the quantum 5" matrix in the same manner as the classical, by taking the 
massless limit {G 0) of the sine-Gordon S matrix [|5|. The Lagrangian for the massive 
theory is 



{df,(pf + Gcos, 



+ g cos 



(3.1) 



The physics of the model depends crucially on the parameter (3, which can be scaled out 
when studying the classical equations of motion. It is convenient to define the parameter 

In the classical limit jS 0, the results of this section of course reduce to those of the 
previous section. We should note that a massless 5" matrix can often be derived directly, 



without taking a limit of a massive theory |T2,13,11.17]. This is done by imposing the 



constraints of the Yang-Baxter equation, crossing and unitarity. 



3.1. The massless bulk S matrix 



The bulk structure of the quantum sine-Gordon theory is very well known P,|2^,p5 



At any value of /3, the spectrum includes a soliton S and an antisoliton A of mass m, with 
m oc These are the quantized particles corresponding to the classical solitons 

discussed in the previous section. Moreover, we have bound states (breathers) indexed by 
n a positive integer less than A, with mass = 2m sin nn /2X. As for the classical theory, 
we obtain the massless limit by sending the rapidities of these particles to ±oo and at 
the same time scaling appropriately the mass m — > so that the energy of each particle 
remains finite. In this limit we obtain a set of left-moving and right-moving particles. The 
right movers have dispersion relation E = P and can be parametrized as 

Es,A = /ue^ Er, = 2^ sin e^ (3.2) 

where fi and 9 are now renormalized mass and rapidity. The left movers have Es^a — 
— P = |Ue~^ and likewise for the breathers. 

The left-left and right-right two-particle S matrices are given by the same formula as 
in the massive case. The matrices Sll and 5"^^^; depend on the difference of renormalized 
rapidities because the only relativistic invariant that can be formed is the ratio of momenta. 
(There are no dimensional parameters in the bulk problem). The S matrix is purely 
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elastic; individual momenta do not change in a collision. Although the classical S matrix is 
diagonal, the quantum one is not: the initial state \S{9i)A{92)) can scatter to \A{9i)S{92)) 
because A and S have the same mass. For the soliton and antisoliton one has the three 
usual amplitudes 

a{9) = sinX{7r + i9)Z{9) 

b{9) = -sin Xi9Z (9) (3.3) 
c{9) = smXnZ{9). 

where the element 0(6*1 — 6*2) describes the process \S{9i)S(92)) \S(9i)S(92)) , 6 describes 
SA SA, c describes the non-diagonal process SA — > AS, and there is a symmetry under 
interchange of soliton to antisoliton (corresponding to ^ ^ —(p)- The function Z{9) is a 
well-known normalization factor, which can be written as 

7,,, 1 f - r dy . 29yX sinh(A - l)y 

Z{9) = — — — - exp [i 7^ sm 



sin A(7r -|- i6') \ J_^2y tv sinhycoshAy 

Notice that when A is an integer, the scattering is diagonal and that a = ±b. 

The left-right scattering comes from the 9 ^ 00 limit of (|3.3| ); it is diagonal and 
rapidity-independent, but this constant phase is different for S'5' ^ 5'S' and SA ^ SA. A 



non-constant Slji induces a flow between bulk critical points ||T^,^|Tl[; we will not discuss 
this here. 

The breather-soliton and breather-breather S matrices are well known . The mass- 
less limit is taken in the same way, so Sll and Srr are the same as the massive S matrix, 
while Slr is a constant. 

3.2. The boundary S matrix for solitons 

To study the effect of the dynamical boundary condition we now consider the scat- 
tering of these massless particles off the boundary. Since the boundary introduces a scale 
Tb (the "boundary temperature") to the problem, the boundary S matrix element for a 
particle with momentum P depends on the ratio P/Tb- Defining the boundary scale 9b 
via Tb = /ue~^^, we see that the S matrix element for a left mover to scatter off the 
boundary (and thus become a right mover) depends on 9 — 9b- 

We can obtain the S matrix from the results of 0. This S matrix has two free 
parameters, rj and 0. Setting ^0 = in the boundary potential as we have done corresponds 
to setting 77 = 0. The remaining parameter O is a function of the boundary scale g. The S 
matrix then depends on a product of functions of a. Act -f G, and Act — O, where a here is 
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the unrenormalized rapidity (which is caUed ^ in ) . To get the scattering of left movers 
off of the boundary, we simply take the limit of very large and negative ol and appropriately 
scale the parameter O, so that a — (0/A) = 9 — 9b remains finite. The effect of this scaling 
is to obtain left-moving particles of energy comparable to the boundary term energy. 

There are two different boundary amplitudes: P{9 — 9b) for the S ^ S and A ^ A 
processes, and Q{9 — 9b) for ^ A and A S. The symmetry (f) — * —cj) ensures that 
the amplitudes are the same when soliton and antisoliton are interchanged. This S matrix 
is not diagonal; notice that the boundary perturbation breaks the U{1) soliton- number 
symmetry, thus allowing soliton to scatter into an antisoliton at the boundary. A.s g ^ oo, 
this symmetry is restored, so at the Dirichlet fixed point, the element Q must vanish. 
Taking the limit of the result of as described, one finds 



P(6/) = e^^/^R{9) 
Q(^e) =ie-^'^/^R{9) 



(3.4) 



where the function R reads 



3*7 



YiiO) 



^^^^ 2cosh(f -zf)ny,(_^) 

y^,, ^ r(| + /A-Ag)r(i + (/ + l)A-Ag) 

r(i-K(/ + i)A-A^)r(| + (/+i)A-A^)- 

A useful integral representation of R is given by 

2cosh(^-zf) ^ y_oc 2y TT sinh2|/coshA|/ 

The constant phase 7 depends on A and turns out to be necessary to satisfy the cross- 
unitarity relation . For right movers scattering off a boundary, the answer is similar but 
differs in some constant phases and signs. 

We can check this answer in several limits. As 9b 00, we obtain the Neumann 
fixed point. Since P{9 — 9b) in this limit, the boundary scattering is completely off- 
diagonal, so that the U{1) soliton-number symmetry is maximally violated. As 9b —00, 
we reach the Dirichlet fixed point, where the scattering is diagonal. The results at both 
limits follow from the analysis of 0. 

We can also check the answer at several values of A. The case A = 1, where the 
perturbing operator has dimension | , has already been solved |jT^ . It was shown explicitly 
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in the corresponding lattice model that the bulk theory is equivalent to two critical Ising 
models, and that the boundary interaction corresponds to a boundary magnetic field in 
one of the two Ising models. We have the same result here. The bulk theory corresponds 
to a free massless Dirac fermion, where the soliton corresponds to the fermion and the 
antisoliton the antifermion. This is easily seen in ( |3.3| ) : the S matrix elements are a = b = 
1, c = 0. We can decompose this into the massless Majorana fermions of two critical Ising 
models by looking at the linear combinations ( | S") + 1 A) ) and {\S) — \A)). The corresponding 
boundary S matrix elements are then P + Q and P — Q. For A = 1, we then have 

P + Q = l 

^ ^ , fO in 
P ~Q = tanh 



2 4 

The S matrix for an Ising model in a boundary magnetic field indeed is tanh (f — as 
can be seen by explicit calculation or by taking the limit of the massive S matrix in . 

Another interesting value is A = 0, where the dimension of the perturbing operator 
is 1, so that the operators exp(i47r(/)), exp(— z47r(^) and dcj) are the generators of an SU(2) 
symmetry at the boundary. The fact that cos47r^ can be rotated into dcf) means that the 
operator is exactly marginal, and the boundary interaction does not break the conformal 
symmetry. The resulting boundary S matrix cannot depend on rapidity: because there are 
no dimensional parameters in the theory, there is no way of making a relativistic invariant 
out of a single momentum. Indeed, we see that the rapidity dependence in the S matrix 
( ^.4| ) goes away. This does not mean that the S matrix is trivial since we can scale 6b with 
A. Actually, in the limit A ^ this is necessary to get a finite value of g, following the 
relation of g and 6b (see (4.16) below). Returning to the original parameter O, we have 

ee/2 

P 



V 2 cosh O /„ „N 

,-e/2 (3.6) 

V2 cosh e 

This is now a simple rotation in the 5", A space 

cos -ngii i sin ngu 
i sin ngn cos Trgn 

with 

tanngji = e~®. 

This is exactly the picture obtained in [p6| , |27| , p8| , where gn is the renormalized coupling 
constant. 



(3.7) 
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3. 3. The boundary S matrix for breathers 

One can easily treat similarly the boundary S matrix for the nth breather The 
integrability requires that it be diagonal. One finds amplitudes that depend on the parity 
of n: 

i i In {e)=2y cosh^cos(/-l/2)./A _ 
%de ^ cosh2 6'-sin2(/- l/2)7r/A 

and 

lAini2(2't-i)m = ^ I 2V cosh 6> cos /tt/A 
ide ^ ' cosh e ^ cosh^ 6 - sin^ I-k /A ' 

The integral representation is 

o(n)/m ■ dy . 2Xey sinhny 

'(6)=expi / — sm — — . (3.8) 

J_^2y n smhj/coshA?/ 

4. The exact boundary free energy 

Using this scattering description we can now compute the evolution of the boundary 
entropy between Neumann and Dirichlet boundary conditions. Recall that if we consider 
the theory defined by the Lagrangian ( |1 . 1| ) at temperature T the free energy reads 



F = fL-Tln{gig2), (4.1) 

where the bulk term / is L independent, Ingi and In (72 are boundary entropies associated 
with the two boundaries at a = 0, L. The boundary interaction does not affect /, but it 
causes gi to fiow from its Neumann value to its Dirichlet value. As we will see, gi decreases, 
in accordance with the "(7-conjecture" of |^^. In this section, we explicitly compute gi as 



a function of the boundary scale by using the thermodynamic Bethe ansatz (TBA) P,pi] 

4.1. The TBA equations in the presence of a boundary 

The TBA exploits the fact in an integrable model, we can find the exact relation 
between the density of states P and the particle density p. (Often the quantity P — p 
is called the density of holes p^ .) In models that are not integrable, one generally uses 
the free-particle density of states and then uses perturbation theory in the interaction 
parameters. However, the fact that the scattering is completely elastic and factorizable 
in an integrable model gives an exact functional relation between P and p, known as the 
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Bethe equation. We then write down the free energy in terms of P and p. The particle 
density p of the system at temperature T is then found by minimizing the free energy, 
using the relation of P to p. 

We discuss the case when the bulk and boundary S matrices are diagonal. When there 
are p different species of particle, we have the two-particle S matrix elements Srs{Oi — 6*2)5 
where r and s run from 1 to p. These S matrix elements are the phase shift in the 
wavefunction when two particles are exchanged. We also have the boundary S matrix 
elements Rr{6 — 6 b) , which gives the phase shift when a particle of species r bounces off 
a wall and changes its rapidity from 6 to —6. We have a gas of M particles on a line of 



length L, with the zth particle of species r^. As in the Kondo problem [|T3, we map the 
problem onto a line of length 2L (— L < a < L) by considering the right movers to be left 
movers with cr < 0. Thus we have only left movers scattering among themselves and off of 
the boundary, which can now be thought of as an impurity (a particle with rapidity Ob)- 
(This trick is common to boundary conformal field theory and can only be used in the 
massless limit.) For simplicity, we put periodic boundary conditions on the system; these 
do not change the boundary effects at a = 0. 

First we write the Bethe equations. These equations quantize the set of allowed 
rapidities for a system of left-moving particles on a circle of length 2L with an impurity 
at cr = 0. They are obtained by collecting all the phase shifts due to particle-particle 
and particle-impurity scattering when arguments in the wave functions are analytically 
continued a ^ a + 2L. Demanding that the wavefunction be periodic gives the constraint 

M 

where E = —P = prexp{—9) for a particle of type r. Since the scattering is diagonal, we 
can define the densities pr and for each species of particle, so that {pr{d) + p^{0))2Ld6 
gives the number of allowed rapidities between 6 and 9 + d9 for species r, and pr gives the 
density of filled states. Taking the derivative of the logarithm of ( [4.2|) gives one equation 
for every type of particle: 



P 1 

27r{pr{e) + p'^M) = l^re-' + J2^rs* Ps^O) + —^^{9 - 63), (4.3) 
where -k denotes convolution: 

/oo 
da'f{a - a')g{a') 
-00 
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and ^ 

Kr{e - Ob) = -i-^lnRrie - Ob). 

uu 

The effect of tlie boundary is seen in tlie last piece of ([4.3D, proportional to 1/L. 
We now consider tlie free energy T = E — TS. For the energy we have 

£ = 2L de^pr{e)nre-^ 

r=l 

r ^ 

S = 2L deJ2 [iPr + Pr) Hpr + Pr) " Pr In - In p^] . 
•J r=l 



(4.4) 



(4.5) 



and entropy 



(4.6) 



(The particles act like fermions when filling levels.) The free energy is the value of —T In Z 
at the saddle point of the partition sum or path integral: as a function of p, JF is at a 
miminum. Via 

dJ^ dT ^ dT ap^ ^ 
= V > — = 

dpr dpr ^ dp^ dpr 

one obtains equations for densities, using (|4.3| ) to determine dp^/dpr- These equations 
do not depend on the boundary term, because boundary terms only appear in the Bethe 
equations, and disappear in a variation of the densities. Defining 

Pr{ey 



one finds the well-known equations PJTO| 



^(^) = ^^e-' - J- V ★ ln(l + 

U ZTT ^ — ^ 

s=l 



(4.7) 



where e(^) = e{6 + ln(p/T)). Boundary terms enter when expressing the free energy in 
terms of the e^. We find them by rewriting the energy ( [4.5| ) as 

8 = 2L Y.Pr 



Ter + ^ <^rs * ln(l + e 

s=l 

. r P P / 

^LT / J2Pr^r + J2(P-+P 

Lr=l s=l ^ 

P 



h _ P±^-B _ J_ ^ 

' 2tt 4ttL V 



ln(l + e-"^) 



r=l 



TLpr -e ,T 

TT ZTT 



ln(l + e" 
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where we used (|4.7| ) to get to the first hne, (|4.3|) to get to the second, and ( [4 .61) to get to 



the third. Thus we have 

:F = T,uik-T I ^^«:,(^-ln(T/TB))ln(l + e-^' (^)). (4.8) 

It is a general result that Tbuik = — ^T'^L in a massless bulk theory, where c is the central 
charge of the conformal field theory ||31|| . 

Although the equations (|4.7| ) for e{9) cannot be solved explicitly for all temperatures, 
the free energy is easy to evaluate as T ^ and T oo, as we will show in sect. 
4.3. Moreover, one can extract the analytic values of critical exponents by looking at the 
form of the expansions around these fixed points. Also, they are straightforward to solve 
numerically for any T. 

Several notes of caution are necessary. At the order we are working, the formula 
( [4.6|) for the entropy is not quite correct, because there are 1/L corrections to the Stirling 
formula used in its derivation. Also, at this order, the logarithm of the partition function 
is not S — TS: it depends not only on the saddle point value of the sum over all states, 
but also on fiuctuations. Their net effect is that we cannot compute the g factors from 
J-' alone. However, both of these corrections are subleading contributions to the bulk free 
energy, and do not depend on the boundary conditions. Therefore we can still compute 
differences of g factors from JF; the corrections are independent of the boundary scale 9b 
and cancel out of the difference IST 



4-2. The exact boundary free energy for X integer 

In this section we consider only the case A a positive integer. At these values, a major 
simplification takes place: the bulk scattering ( p.3|) is diagonal. The boundary scattering 
still is not, but we can redefine our states to be |±) = {\S) ± \A))/V2 as in the previous 
section, so that the boundary scattering is now diagonal: 



R+{e) = P + Q = cosh 

R-(^0) = P_Q = cosh (^^ + e-'^/^R{9) 



xe 


in 


Y ~ 


~4 


xe 


in 



^ If A is even, this actually makes the bulk scattering completely off-diagonal (e.g. | + +) 
scatters to | )), but the TBA equations turn out the same. 
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Using the analysis of the previous subsection, we can read off the answer. There are 
A + 1 particles: the breathers, which we index n = 1 ... A — 1, and the ± particles. We 
define the Fourier transform as 

fiy) = r ^e^''''^-m. 

J-oo 

The bulk kernels (p are well known [^^; they can be written in the form 



coshycosh(A - n)ysinh/?/ 

fni = Oni - 2 — — — n, / = 1 . . . A - 1; n>l 

cosh Ay smh y 

cosh y sinh ny 

^' cosh Ay sinh y 

sinh(A — l)y 



2 cosh Xy sinh y 

with (frs = fsr- The boundary kernels follow from ( [4.4| ), and are 

sinh ny 



2 sinh y cosh Xy 



sinhfA — l)y 1 , . . 

K- = H r- 4-10 

2 sinh 2y cosh Xy 2 cosh y 

sinhfA — l)y 



2 sinh 2y cosh Xy 

The equations for the e^. can be written in a much simpler form by using a few trigonometric 
indentities. One finds 

er = Ki.^NrsMl + e^^), (4.11) 

s 

where the label s runs over breathers and ±, the kernel K{6) = l/27rcosh^, and Nrs is 
the incidence matrix of the following diagram 

12 s A-2 / 

O O 0---0 O A-1 



The dependences on the mass ratios seems to have disappeared from (|4.11D, but they 



appear as an asymptotic condition: the original equations (|4.7|) indicates that the solution 
must satisfy 

er — e as y — ^ — oo. 
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To recapitulate: the impurity contribution to the free energy is 

T.mp = / 1^ E ^-(^ - M^/Ts)) ln(l + e-^'-(^))), (4.12) 



where the are given by (|4.10| ), and the e^. are the solutions of the non- linear integral 
equations (|4.11| ). 



At non-integer values of A, the bulk S matrix is not diagonal. As a result both bulk 
and boundary scattering have to be simultaneously diagonalized. It is known how to do 



this in the bulk |3^], and the inclusion of the boundary can be presumably be acomplished 
by using the quantum inverse scattering method along the lines of |^5[ and [^2|. The 
simplest values to do should be A = 1/t with t integer, where there are no breathers, and 
the bulk TBA system is given by ( [4.11| ) with different conditions for ^ — > — oo. In fact, 
using the ideas of [^2| and the "duality" of TBA equations for X = t and A = y we have 



been able to conjecture a TBA for the latter case that reproduces all expected results for 
our flow. For simplicity we do not present this TBA here. 

4-3. The Neumann and Dirichlet limits 

We can evaluate the impurity free energy explicitly in several limits. In the IR limit 
T/Tb — > the integral is dominated by — — oo where the source terms in (f4.7| ) become 
very big. Hence er{—oo) = oo and the impurity free energy vanishes in this limit. In the 
UV limit T/Tb — > oo the integrals are dominated by the region where 9 is large so that 



the source terms disappear in (|4.7| ) and the go to constants: 



= e'"(°°) = (n + 1)2 - 1; x± - A. (4.13) 



Therefore we obtain 

9N imp 



In 



9d T 

t-2 



imp 



T 

UV 



IR 



(4.14) 



J2 ^^"^ ln(l + 1/xr,) + (/(+) + /(-)) ln(l + l/x±), 



n=l 



where 



^^^^^ I ^^r{e) = m- 
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Thus = n/2 and + = A/2, and we find 



ln^ = ^ln^ + y ^In , , 

gn 2 A ^^4 n(n + 2) ^^.^^^ 

= ^ln(A + l). 

This is agreement with the ratio calculated from conformal field theory, as detailed in the 
Appendix. 

We can also find the dimension of the perturbing operators. From the equations ( [4.7| ) 



one deduces the following expansions for T/Tb large: 

3 

As a result it is straightforward to see that near (7=0, JF can be expanded in powers of 
(Tb/T)^^/*^^+^). On the other hand we expect JF to be an analytic function of g"^. Hence 

go,{y.e-'-)^'^^^'\ (4.16) 

This agrees with the conformal result that the perturbing operator cos[/30(O)/2] has bound- 
ary dimension a; = 1/(A + 1) = /3^/87r. In the IR limit of T/Tg small, one can expand 
out the kernels Kj- in powers of exp(^;, — 6). This leads to the fact that the irrelevant 
operator which perturbs the Dirichlet boundary conditions has dimension x=2. This is 
the energy-momentum tensor. (We note that there is another irrelevant operator in the 
spectrum with dimension x = A 4- 1, which for < A < 1 is the appropriate perturbing 
operator [||.) 



5. Future directions 

We see that at A = 0, the boundary 5" matrix is a constant for all 9. This is expected, 
since we can continuously interpolate from Neumann to Dirichlet without destroying the 
conformal symmetry. Using the conformal symmetry, one in fact is able to calculate the full 
partition function along this interpolation |]26| - |2^] , not just at the Neumann and Dirichlet 



points as in the Appendix. It would be interesting to precisely understand how to map 



the scattering basis discussed in to the scattering basis used here. This would provide 



a great deal of insight into the long-confusing relation between the states arising from the 
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Bethe Ansatz to those in the Virasoro towers of a conformal field theory (see |3^, for 
example) . 

The quantum 5' matrix for Dirichlet boundary conditions ^(0) = can be derived 
explicitly by applying the Bethe ansatz to an underlying lattice model, the XXZ spin chain 
with a boundary magnetic field [0. Presumably the more general 5" matrix studied here 
corresponds to the XXZ chain with boundary terms of the type that break the U{1) 
symmetry [EH . 



By expanding the partition function of our system in powers of g we obtain a one- 
dimensional neutral Coulomb gas of particles on a circle, whose solution is therefore indi- 
rectly provided by the TBA. Previously, such a Coulomb gas had been solved for only one 
type of charge. 

Perhaps the most interesting future direction is the application of these results to the 



experimentally realizable system of tunnelling between fractional quantum-Hall edges ||3^ . 



We will report on this application soon |4C 
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Appendix A. Neumann and Dirichlet partition functions 

We here use standard conformal theory techniques to calculate the partition function 
of our quantum system on a cylinder of length L and circumference 1/T and Neumann or 
Dirichlet boundary conditions at the ends. 

With Dirichlet boundary conditions on both sides, winding modes are not allowed 
around the cylinder, but they are allowed along it. Thus the field (j) obeys ^(a = 0) = 
cj){a = L) -\- n4n / f3 with n an integer. (In the usual conventions, this corresponds to a 



bosonic radius of v47r//3.) Hence |^ 



1 V^8.nV/3^ (A.l) 



where w = e '^/^'^ , After modular transformation one finds 

1/2 
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where q = e"^'^^^. With Neumann on one side and Dirichlet on the other, no wind- 
ing modes of any sort are allowed. The energy, however, is quantized with half-integer 
eigenvalues. Using the Jacobi triple product formula gives 



ZND ^^y^in-l 



/2f/4 



(A.3) 



which after modular transformation is 



V2r]{q 



(A.4) 



By taking the ratio in the limit of large L, we find 
in agreement with (|4.15|) . 



(A.5) 
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